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H 要 : 随 着 油气 资源 需求 不 断 增加 ,水 平井 应 用 越 来 越 广 泛 ,水 平井 钻 柱 的 横向 振动 将 直接 影响 
钻井 效率 和 钴 柱 安全 。 与 直 间 不 同 ,水 平 并 存在 造 斜 段 和 水 平 段 , 钻 柱 的 初始 变形 和 受 力 条 件 更 加 
复杂 。 为 了 准确 分 析 水 平 并 钻 柱 的 横向 振动 ,文章 综合 考虑 三 维 并 眼 轨迹 、 钻 并 参数 和 钻 柱 与 并 辟 
接触 等 参数 ,基于 Lagrange 方程 建立 了 水 平 并 钻 柱 非 线 性 动力 学 特性 数学 模型 ,采用 有 限 单元 法 对 
之 模型 离散 并 求解 ,并 选取 四 川 某 水 平 并 实际 参数 进行 验证 ,最 后 讨论 了 钻 压 、 转 速 和 扶正 器 对 水 平 
回 并 钻 柱 横向 振动 的 影响 。 研 究 结果 表明 :在 钻 进 过 程 中 , 直 并 段 钻 柱 横向 振动 较 小 , 造 斜 段 和 水 平 
人 7 段 钻 柱 横向 振动 较 大 , 越 靠近 钻头 处 ,横向 振动 越 大 ; 钻 压 对 水 平 段 和 造 儿 段 钻 柱 横向 振动 有 较 大 
GO 影响 , 钻 压 越 大 ,该 段 钻 柱 横向 振动 越 大 ,而 直 并 段 钻 柱 受 钻 压 影响 很 小 ;在 文章 分 析 条 件 下 (转速 
€ 为 30 ~ 120 r/min) ,全 并 段 钻 柱 横向 振动 无 明显 变化 ;扶正 器 可 以 限制 水 平 段 钻 柱 的 横向 位 移 , 且 
CN 扶正 器 安装 位 置 对 其 有 很 大 影响 。 文 章 所 提出 的 模型 也 适用 于 复杂 结构 并 钻 柱 动力 学 特性 研究 ， 
二 研究 成 果 为 水 平 并 钻 并 提 效 和 钻 柱 安 全 提供 理论 指导 。 
-关键 词 : 钻 柱 动力 学 ;横向 振动 ;水 平 并 ; 并 眼 轨迹 
ip 中 图 分 类 号 .TE24 文献 标志 码 :A DOI:10. 11776/j. issn. 1000-4939. 2023.03. 022 


Analysis on the lateral vibration of drillstring in horizontal wells 
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Abstract : With the increasing demand for oil and gas , horizontal wells are used more and more frequently. 
The lateral vibration of the drillstring in horizontal wells will directly affect drilling efficiency and safety. 
Different from vertical wells, horizontal wells have buildup and horizontal sections besides vertical section; 
thus the initial deformation and stress conditions of the drillstring are more complicated. In order to accu- 
rately predict the lateral vibration of the drillstring in horizontal wells , based on the lagrange equation ,a fi- 


nite element model of drillstring system is established with consideration of the 3D borehole trajectory, 
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drilling parameters and contact between the borehole and drillstring. The model is discretized and solved 
by the finite element method. Then, the lateral vibration characteristics of the drillstring are verified by the 
field data of a horizontal well in Sichuan. Finally ,the influence of the weight on bit , rotation speed and sta- 
bilizer is analyzed. The results show that in the drilling process , the lateral vibration of the drillstring in the 
vertical section is slight; the drillstring lateral vibration in the buildup and the horizontal sections is vio- 
lent, and the closer to the bit is, the more violent the lateral vibration is; the weight on the bit has small 
effect on the lateral vibration in the vertical section , but it has significant effect in the buildup and horizon- 
tal sections. the larger the weight on bit,the more violent the lateral vibration of the drillstring. Under the 
actual drilling conditions (the rotation speed is 30 — 120 r/min) ,the change of the lateral vibration of all 
drillstring sections are not obvious. The stabilizer can limit the lateral displacement of the drillstring in the 
horizontal section , and its location can greatly affect the stabilizer effect. The model can also be used to 


study the dynamic characteristics of drillstring in cured wells , and the results can provide theoretical guid- 


ance for improving drilling efficiency and ensuring drillstring safety. 


.. Key words : drill string dynamics; lateral vibration ; horizontal well;3D borehole trajectory 


目前 ,为 了 开发 页 岩 气 藏 等 特殊 油气 藏 和 提高 
酒 绊 单间 产量 ,水 平井 应 用 越 来 越 广泛 "”) 。 水 平井 
FIERE , 钻 柱 在 窗 ,长 井 眼 内 受到 轴 向 力 ,横向 力 和 
摩 溺 力 等 载荷 联合 作用 ,会 产生 轴 向 振动 .横向 振动 
RUDI a A SEC 。 其 中 ,横向 振动 被 认为 
是 危害 最 大 的 振动 类 型 5 :@ 钻 柱 对 井 壁 的 高 频 冲 
起 异常 严重 , 井 壁 失 稳 、 掉 块 卡 钻 和 落石 垮塌 等 井下 
KRRB O ;@) 易 造成 钻 柱 磨损 ,裂纹 , 刺 漏 变 
形 等 ,使 得 钻 柱 失效 !9 ;@) 钻 头 侧 向 力 和 钻头 转角 在 
钼 特 与 井 壁 碰撞 时 会 发 生 较 大 变化 ,导致 实际 井 眼 
锦 迹 发 生 偏 移 ") 。 这 些 危害 在 水 平井 钻井 中 尤为 突 
出 5 但 是 ,在 水 平井 钻 进 时 ,底部 钻 柱 向 井口 传递 时 
机 丙 振 动 有 很 大 减弱 ,即使 井 底 钻 柱 横向 振动 剧烈 ， 
也 全 法 在 地 面 上 显示 出 来 ”。 因 此 ,水 平井 钻 柱 的 
横向 振动 研究 ,对 防止 钻 柱 失效 .提升 钻井 效率 和 控 
制 井 眼 轨迹 有 重大 意义 。 

针对 钻 柱 动力 学 行为 ,国内 外 众多 机 构 和 学 者 
已 经 开展 了 大 量 分 析 。 其 中 ,高 德 利 等 ”通过 构 
建 三 维 非 线性 动态 控制 方程 来 分 析 钻头 处 的 钻 具 组 
合 动力 学 行为 ,并 研究 了 压 扭 作 用 下 该 处 钻 柱 的 届 
曲 特性 。 刘 清 友 等 "中 基于 弹性 杆 理 论 建 立 的 钻 柱 
动力 学 模型 ,包括 了 岩石 和 牙 轮 钻 头 之 间 的 作用 。 
李子 丰 等 5 系统 研究 了 钻 柱 的 稳定 性 问题 ,包括 
斜 直 井 中 钻 柱 的 几何 线性 、 非 线性 螺旋 届 曲 和 稳定 
性 判别 。 管 志 川 等 ”为 了 分 析 底部 钻 具 的 动力 学 
特性 ,制作 了 BHA 模拟 装置 ,并 对 钻 柱 井 底 的 运动 
情况 和 振动 波 传播 进行 了 分 析 。 关 铁 等 ”2 选择 斜 
直 井 眼 整体 销 柱 为 对 象 ,分 析 了 钻 柱 的 轴 向 振动 ,并 


计算 其 固有 频率 。 况 雨 春 等 利用 ANSYS 软件 ， 
针对 不 同 底部 钻 具 建立 了 特征 值 届 曲 分 析 方法 ,得 
到 其 临界 失 稳 载荷 ,为 防 斜 打 快 提供 理论 基础 。 祝 
效 华 等 构建 的 直 井 钻 柱 静 力学 模型 和 其 耦合 振 
动 模型 ,用 于 分 析 钻 井 液 钻 井 和 空气 钻井 时 钻 柱 振 
动 和 受 力 。 狄 勤 丰 等 站 基于 转子 运动 与 受 力 关系 
来 研究 下 部 钻 具 组 合 的 动力 学 行为 ,通过 数值 仿真 
研究 其 涡 动 特性 。 付 建 红 等 ”为 了 解决 钻 进 过 程 
中 井 径 扩大 的 问题 ,基于 弹性 力学 原理 建立 了 侧 钻 
井 钻 柱 动力 学 模型 。 练 章 华 等 的 根据 气井 水 平井 
钻 柱 的 实际 结构 和 相似 原理 ,研制 相关 实验 装置 , 通 
过 直 梁 有 限 元 模型 分 析 了 200 m 水 平井 底部 钻 具 的 
届 曲 行为 和 接触 作用 。 刘 强 等 ”1 利用 动力 学 软件 
Adams 对 螺杆 钻 具 的 钻 柱 动力 学 特性 进行 了 研究 。 
邓 继 华 等 在 ANSYS 软件 中 ,开发 了 一 种 几何 非 
线性 平面 梁 单 元 ,有效 缩短 了 计算 时 间 。 张 毅 成 
等 5 利用 连续 体 振动 理论 ,开发 了 多 维 激励 下 的 梁 
结构 响应 求解 方法 。 杨 超 等 5 通过 算 例 ,对 比分 析 
了 5 种 积分 算法 对 非 线性 振动 系统 的 响应 情况 。 杜 
文静 等 2 通过 ~ 级 多 步 块 方法 对 Hamilton. 系统 进 
行 了 求解 。 另 外 ,国外 学 者 也 对 此 开展 了 大 量 研究 。 
HANSON 等 考虑 钻头 与 井 底 的 相互 作用 以 及 钻 
头 能 量 损 失 等 ,建立 了 井 底 钻 柱 的 动力 学 模型 ,来 研 
究 PDC 钻头 的 性 能 。KHULIEF 45?! 基于 Lagrange 
方程 建立 了 钻 柱 的 动态 方程 组 ,通过 变换 得 到 其 降 
阶 模型 ,并 应 用 该 模型 计算 了 1200 m 垂直 钻 柱 的 模 
型 特征 。GULYAYEY 等 ”考虑 流体 阻力 扭矩 、 轴 
向 力 等 对 钻 杆 动态 特性 的 影响 ,研究 了 钻 柱 的 固有 
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频率 和 振 型 。ZARE 等 55] 借助 对 钻 柱 的 非 线性 静 力 
学 研究 ,确定 其 有 限 振动 长 度 。BESSELINK 等 为 
了 探索 导致 务 滑 振动 的 机 理 , 对 此 建立 了 钻 柱 轴 向 - 
扭转 振动 的 动力 学 模型 。STROUD 457" 开发 了 一 
个 具有 2 个 自由 度 的 解析 模型 ,用 于 预测 直 井 钻 具 
组 合 的 动态 特性 。VLAJIC 450 等 构建 了 一 个 4 自 
由 度 的 降 阶 方程 ,来 获得 垂直 钻 柱 在 旋转 运动 期 间 
的 扭转 振动 特征 。NGUYEN 457 使 用 粱 有 限 单元 
法 对 钻 柱 动力 学 进行 了 时 域 建 模 和 验证 ,该 模型 可 
以 考虑 3 000 m 外 有 具 组 合 的 真实 几何 形状 。SARK- 
ER 45" 考虑 钻 柱 与 井 壁 接触 时 的 黏 滑 和 旋转 作 
用 ,建立 了 非 线性 三 维 多 轴 模型 ,用 于 预测 水 平井 底 
部 钻 具 的 振动 情况 。TRAN 4& ^ aep T ski ts god 
的 相互 作用 和 钻 柱 与 井 壁 的 接触 作用 ,采用 Timosh- 
enko 梁 单 元 法 建 模 , 分 析 了 钻井 液 的 影响 。 

Cr E, E RUIRALAIMEORETEE ,水平 并 和 定向 
指引 柱 的 动力 学 行为 进行 了 大 量 研究 ,本 研究 在 此 
基础 上 主要 分 析 了 钻井 参数 对 水 平井 钻 柱 横向 振动 
ita. 本 研究 综合 考虑 三 维 井 眼 轨迹 、 钻 具 组 合 、 
tibus 等 参数 ,基于 Lagrange 方程 建立 了 水 平 
AER HEAR 性 动力 学 特性 数学 模型 ,通过 有 限 
部 克 法 离散 并 求解 ,最 后 选取 四 川 某 水 平井 实际 并 
Wess 与 钻井 数据 ,对 其 非 线 性 横向 振动 动力 学 行 
X T WE. 


kk 平井 钻 柱 系统 动力 学 模型 的 建立 


二 水 平井 由 3 部 分 组 成 ,包括 直 井 段 , 造 斜 段 和 水 
平 蜂 ( 图 1) 。 直 井 段 钻 柱 上 部 主要 受 转 盘 扭矩 和 大 
钩 载荷 , 造 斜 段 和 水 平 段 钻 柱 由 于 重力 影响 ,频繁 与 
井 眼 下 端 接触 , 钻 进 时 除了 受 井 壁 间 摩擦 力 影响 ,还 
会 受到 泥浆 阻尼 的 作用 。 


XE TT 


图 1 水 平 并 
Fig.1 Horizontal well 
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因此 ,本 研究 在 分 析 钻 柱 动力 学 行为 时 ,做 出 以 
下 假设 :中 视 井 简 为 等 截面 圆 , 井 眼 轨迹 连续 光滑 ; 
@ 将 销 柱 视 为 具有 均匀 材料 和 几何 特性 的 三 维 弹 性 
梁 , 其 变形 在 弹性 范围 内 ;(3) 忽 略 钻 柱 之 间 的 连接 螺 
纹 和 接头 ;(4) 将 扶正 絮 视 为 大 尺寸 钻 柱 。 

采用 双 节 点 梁 单 元 对 水 平井 眼 连 续 钻 柱 进行 离 
散 ( 图 2) 。 每 个 单元 有 6 个 自由 度 ,x、y、z 表示 平移 


量 ,9、6. 表示 横向 旋转 角 ,6, 表示 扭转 角 。 钻 柱 运 
动 可 表示 为 
U, = [ xi ,yi ,2i1, 0 ,0,04,%;,);, 


2; ,8,,0,,0,] 


EK2 节点 单元 
Fig.2 Element node 
控制 钻 柱 运动 的 Lagrange 方程 可 表示 为 
[中 人 LT e Y) 
dt 9 U, ðU; 
式 中 :7 表示 动能 ;V 表示 势能 ;上 面 的 点 表示 变量 
U, 对 时 间 的 导数 ;F, 表示 施加 在 节点 上 的 外 力 。 
梁 单 元 的 总 动能 表达 式 为 


Lg 
T = el PAL Cu, 0? + Cu, D? + Cu? 1dx + 
" 1 


=F, (2) 


le 
| pAel2( ,cos( — ù sin(QD dx + oleg 23 
; 


1 [^ l . TM 
3] ^u. * 20, 6,00 * L.C; +0)1dx (3) 


REP: iy o iy à, HS XY Z 轴 上 的 平移 速度 ， 
si 0, x O, x 0 IDRE RENTER TTE x Hy 轴 和 
ARA ond à b, 分别 表示 梁 单 元 截面 线 
re m RR. 
m; 是 指 井下 钻 柱 的 实际 转速 ,rad/s;p 是 钻 柱 的 密 
度 ,kg/m ;1 是 梁 单元 的 长 度 ,m;4 是 梁 单元 的 横 截 
面积 ,mr ;1 是 粱 单元 截面 的 极 惯性 矩 ,mt4; 1 是 染 
单元 横 截 面 的 惯性 矩 ,me; ERG) 中 ,前 3 项 表示 
其 平移 动能 ,后 2 项 表示 其 转动 动能 。 
梁 单 元 的 总 势能 表达 式 为 


投稿 网 站 :http://cjam. xjtu. edu. en 微 信 公 众 号 :应 用 力学 学 报 


- EA 


Mir, 


le - GI, 
(=) dx + — (2) a 
o \ ox 2 ðx 


EIL, ç" 30. ay , Eef 00 
zl.) 71, d 


- PE 
NES + f S t8) + (,*]à dv + 


2 2 
HAE) (2) Je 
w 

2 z^ (6, 2 - 6. Pda (4) 


式 中 :五 是 钻 柱 的 弹性 模 量 ,Pa; G 是 钻 柱 的 剪 切 模 
量 ,Pa。 在 式 (4) 中 ,线性 刚度 矩阵 可 以 从 前 4 项 推 
导出 ,与 钻 柱 弯曲 形变 和 轴 向 形变 相对 应 的 非 线 性 
高 和 矩阵 可 由 第 5 ~ 9 项 推导 出 ,至 于 钻 柱 扭转 变形 
RIS Ill dE 形 对 应 的 非 线性 刚度 矩阵 可 由 最 后 3 项 推 
Sek. 
@ 梁 单元 节点 会 受到 2 种 外 力 ,分 别 为 离心 力 和 
£2. 其 中 梁 单 元 在 x、y、z 上 的 重力 分 量 可 表示 为 


qx = qcosa 


q, = qsina (5) 
q, =0 
su. 4 表示 1 m 钻 柱 的 等 效 重力 ,Nm; a 是 梁 单 元 
轴线 和 垂直 方向 的 夹 角 。 因此 ,重力 矢量 的 等 效 节 
Agios 


《 -gL-gL -gL 
é “F D = qd. dy qy 
EP [=== o o 0 2 
üa q,L ES q,L gL 
-— e NE 全 | (6) 


在 钻 柱 横 截 面 上 ,重心 和 形 心 并 不 完全 相同 , 离 
心力 在 钻 柱 旋 转 时 存在 。 梁 单元 在 x、y、z 上 的 自转 
产生 的 离心 力 可 表示 为 

f, = qel? cos 
i (7) 
jf.=0 
RP p 是 重心 的 相位 角 ,rad。 同 理 ,对 于 节点 的 离 
心力 矢量 ,其 等 效力 可 表示 为 
SHORE: 


2 2 12 
-fL -fL — fÈ -fr 

i 2 2 i 2 =] 
(8) 
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将 式 (3) (4) (6) (8) 代 入 式 (2) ,得 到 
M, U, +C, U, + K,U, = F, (9) 

RP: U, ,U, , U, | F, 表示 随 体 坐 标 系 下 的 广义 加 
速度 .速度 ,位移 和 外 力 矢 量 ; M, LC, K, 分别 表示 
随 体 坐标 系 下 的 质量 矩阵 .阻尼 矩阵 和 刚度 矩阵 。 

式 (9) 为 各 单元 随 体 坐标 系 下 的 运动 方程 ,但 
在 水 平井 实际 钻 进 过 程 中 , 井 眼 轨迹 是 一 条 弯曲 的 
空间 曲线 。 因 此 ,在 钻 柱 动力 学 模型 求解 前 ,需要 对 
真实 井 眼 轨迹 和 钻 柱 初 始 构 形 进行 准确 描述 。 然 后 
建立 一 个 大 地 坐标 系 ,通过 坐标 转化 得 到 各 单元 节 
点 大 地 坐标 系 下 的 运动 方程 。 最 后 着 加 各 梁 单 元 大 
地 坐标 系 下 的 运动 方程 , 即 为 整个 钻 柱 的 运动 方程 。 
1.1 井 眼 轨迹 与 水 平井 钻 柱 的 初始 构 形 

水 平井 井 眼 轨迹 的 基本 数据 可 以 实际 测量 , 利 
用 随 钻 测量 工具 (MWD ) 得 到 。 井 眼 轨迹 坐标 可 以 
根据 方位 角 和 并 斜 角 经 曲率 半径 法 计算 得 到 '* ,最 
后 采用 插值 法 建立 真实 井 眼 轨迹 。 
AX = Rs(sing, - sing) 


AY = Rs Csing, — sing,) 


(10) 
AZ = Rj(sing, - sing, 
AS = RyCsina, - sino; 
As 
R; = 一 一 
/ ^a 
R„( cosa; — cosa, ) 
R; = A 
d (11) 
Aa = oa - Q; 
Ap = p2 - Pı 
As = $s, — 5, 


式 中 : ao 分别 是 井 斜 角 和 方位 角 ,rad;s 是 水 平 投 
影 长度 ,m; Ry 为 井 眼 轨迹 垂直 投影 的 曲率 半径 ,m; 
R; 为 井 眼 轨迹 水 平 投影 的 曲率 半径 ,m; 下 标 “1” 
“2” 为 相 邻 2 个 测 点 。 

当 Aa 和 Ao 存在 等 时 , 则 无 法 计算 曲率 半径 RR， 
和 Rs ,可 通过 以 下 4 种 方法 处 理 。 

a) 当 Aa 20, Ap s OT] UE 
AX = Rs(sing, - sing? 


AY = R,Csing, - sing;? 
AZ = Ascosal 
AS = Assino, 


(12) 
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Assina, _ Assina, 
Ap Ap ^ 
b)4 Aa #0 , Ap =0 时 ,有 


AX = Ryl cosa; — cosa, ) coso; 


其 中 R; = 


= Ry( cosa, — cosa, ) cosQ; 
AY = Ryl coso, — coso, sing, 
1 2 | (13) 
= RyCcosg, — coso, sing; 
AZ = RyCsina, - sina? 


AS = RyCcosa, — cosa; 


其 中 Pr =E, 
c) 当 Aa 20, Ap =0 时 , 井 眼 轨迹 为 直线 ,有 
S AX = Assina, coso; 
m AY = Ascosa,sing, (14) 
e AZ = Ascosa, 
c AB = Aina 


001) * Aa #0 ,{E o, + a, =0 时 ,与 b) 中 情况 相 
lb RU Jr e AB ER RT 

QJ 得 到 测 点 间 的 坐标 变化 量 后 ,通过 坐标 参数 计 
SE Re 每 个 测 点 的 实际 坐标 , 即 


d X, = X, + AX 
s Y, = Y, + AY 
>< (15) 
eg pa eAZ 
E S, = S, 4AS 


1.99 坐标 变换 
坐标 系 转换 关系 如 下 。 


X Cosyz cosyxsinyz  siny,siny; |x 
[ = | — sinyy Cosyxcosyz  siny,cosyz d 
Z 0 — sinyx COSy y Z 


(16) 
式 中 ,yx 和 yz 分 别 表 示 绕 不 和 2 轴 旋 转 的 角度 ， 
rad。 设 T 为 转换 和 矩阵, 即 


cosyz cosyxsinyz sinyxsinyz 
T=| -simy; cosyxcosyz sinyycosyz |(17) 
0 — siny y COSY x 
根据 该 矩阵 ,就 能 得 到 各 单元 在 大 地 坐标 系 的 


运动 方程 ,最 后 通过 合 加 得 到 整个 钻 柱 的 运动 方 


程 , 即 

MU+CU +KU' =F' (18) 
式 中 : Ü U U F 分别 表 示 全 井 钻 柱 上 各 节点 的 
广义 加 速度 矩阵 .速度 和 矩阵、 位 移 和 矩阵 和 外 力矩 阵 ; 
M' ,C' .K' 则 分 别 表示 全 井 钻 柱 的 总 质量 矩阵 .总 阻 
尼 和 矩阵 和 总 刚度 矩阵 。 


1.3 边界 条 件 


1.3.1 井 壁 约束 
在 实际 钻井 过 程 中 , 钻 柱 运动 受到 转盘 的 约束 ， 
其 上 端 贸 接 在 井 眼 中 心 , 承 受 招 矩 和 拉力 。 而 钻头 
T BOURSE HA BCHDERUSU 71 H BB IR DLE S 
到 约束 。 此 外 , 钻 柱 和 井 壁 间 际 较 小 , 钻 柱 径 向 位 移 
会 受 约束 。 图 3 为 钼 柱 和 井 眼 的 接触 关系 。 钴 柱 与 
井 眼 发 生 碰 撞 时 ,将 受到 正 向 力 、 切 向 上 的 摩擦 力 和 

额外 的 摩 氛 扭 矩 。 


(NE A 


5 
5 dj2 o 


H3 ESIR 
Fig.3 Contact between drillstring and borehole 
当 井 简 的 弹性 形变 恢复 时 , 钻 柱 会 反弹 回来 ,而 
节点 又 恢复 自由 。 正 向 力 六 可 以 表示 为 


人 - In( AM 


- [u, - (d, - d;)/2]k, -v, 


, 


T 


ES 


N=) aw >(d,-d;)/2 
0,0 u, « (d, - d,)/2 
(19) 
式 中 :qd, 是 井 简直 径 ,m;d; 是 钻 柱 直径 ,m;v, 是 钻 柱 的 径 
问 速 度 ,m/s;w 是 钻 柱 的 径 向 位 移 ,m;% 是 井 简 的 刚 
度 ,N/m;v1 和 分 别 表示 节点 碰撞 前 后 的 速度 ,m/s。 
切 向 摩擦 力 和 摩擦 扭矩 7 可 以 表示 为 


d, 
F, z BG Fy , T, x FH GO Fy (20) 
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sth, Go.) 是 钻 柱 与 井 眼 接触 时 的 摩擦 系数 ,可 以 
通过 静态 -动力 学 指数 衰减 模型 计算 得 到 , 即 
Mv) = p, + Cu, -pe (21) 
式 中 :uw 和 ,分别 为 动 摩擦 系数 和 静摩擦 系数 ;d, 表 
示 训 减 系数 iu, 表示 滑 移 率 。 
1.3.2 钻头 与 地 层 相 互 作用 模型 
由 于 钻头 与 井 底 岩 石 碰撞 和 其 自身 特性 , 钻 压 
F a 随 着 时 间 + 的 变化 通常 被 表示 为 
Fa) = W + Wsin( wt) (22) 
式 中 : W 表示 钻井 过 程 中 施加 的 钻 压 值 ,N; WW 为 动 
态 钻 压 值 ,N; wj 为 钻 压 波动 系数 ,该 值 同 钻头 类 型 
和 和 钻 柱 转速 有 关 , 即 
o, = njQ (23) 
式 中 : n, 为 钻头 系数 ,对 于 PDC dise n, 21 ,对 于 牙 
AES n, = 3 ; 0 为 钻 柱 转速 。 
DD) 钻头 所 受 摩擦 扭矩 Tu, 表示 为 
Tu) = F Daal Wo + Wsin(w)) (24) 


Dy, 为 钻头 直径 ,m; u 为 地 层 摩 擦 系数 。 


忆 从 上 述 条 件 可 以 看 出 , 钻 柱 动力 学 模型 分 析 较 
为 复杂 ,本 研究 通过 Newmark48 法 ,考虑 碰撞 和 摩擦 
对 模型 求解 ,主要 流程 如 图 4 所 示 。 


输入 基本 参数 : 井 眼 轨迹 、 底 
钻 具 组 合 、 钻 柱 结构 和 钻井 参 


Hp 


计算 等 效 载 荷 矢量 


单元 刚 
阵 和 阻 


Newmark-B 法 计算 外 柱 节点 
位 移 、 速 度 和 加 速度 


图 4 ”有限 元 模型 流程 求解 图 


Fig.4 Flowchart of finite element model solution 
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3 钻 柱 横向 振动 动力 学 分 析 


本 研究 选取 了 四 川 茶水 平井 Al 井 进行 钻 柱 动 
力学 行为 研究 , 表 1 给 出 了 该 井 主要 计算 参数 。 
表 1 Al 井 基 本 参数 


Tab.1 Basic parameters of well A1 


井 眼 长 度 /m 4 200 最 大 倾斜 角 /(? ) 86.2 
钻井 液 密 度 /(kg - m^?) 1200 钻 压 /kN 40 
转速 /(r: min-!) 60 钻头 直径 /mm 165 
钻 狼 外 径 /mm 121 fli EDI fS/ mm 57 
Al EK HE/m 9 [B EEREBE/(kg- m^?) — 7853 
加 重 钻 杆 外 径 /mm 101 加 重 钻 杆 内 径 /mm 65 
加 重 钻 杆 长 度 /m 307 钻 杆 外 径 /mm 101 
钻 杆 内 径 /mm 82 钻 杆 长 度 /m 3 384 
钻 杆 密度 /(kg . m^?) 8676 扶正 器 外 径 /mm 162.5 
扶正 器 长 度 /m 2 


该 并 在 钻 柱 动 力学 模拟 中 , 销 柱 单元 长 度 为 
3 m ,模拟 总 时 间 为 5 s, 时 间 步 长 为 0.001 s, 钻 柱 转动 
5 fl. Fd 5 为 该 水 平井 井 身 结 构图 和 钻 柱 横向 加 速度 
图 。 钴 柱 实际 横向 加 速度 可 以 通过 在 随 钻 测 量 仪 上 
使 用 一 种 有 线 销 杆 技术 ( 距 销 头 36 m) 测 量 得 到 。 


273.05 mm X 1 600 mm 
381.0 mm X 1 600 mm 


196.85 mm X3 600 mm 
241.3 mm X3 600 mm 


127.0 mm X4 200 mm 
165.1 mm X4 200 mm 


(a) 井 身 结构 


CN 


横向 加 速度 /g 
二 一 


时 间 /s 
(b) 横向 加 速度 


图 5 井 身 结构 图 和 横向 加 速度 图 


Fig.5 Wellbore configuration and lateral acceleration 
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从 图 5(b) 中 可 以 看 出 ,该 动力 学 模型 计算 出 的 
钳 柱 横向 加 速度 与 实测 数据 处 于 同一 个 量 级 ,曲线 
基本 吻合 ,因此 该 模型 适用 于 水 平井 中 钻 柱 横向 振 
动 特性 预测 。 

图 6 是 Al 井 钻 柱 轴 向 力 和 横向 变形 图 。 从 图 
6(a) 可 见 ,该 井中 性 点 位 置 大 约 位 于 距 钻头 400 m 
处 。 直 并 段 外 柱 轴 向 力 与 深度 呈 线 性 关系 ,而 造作 
段 和 水 平 段 钻 柱 轴 向 力 的 变化 趋势 随 井 深 增 加 逐渐 
变 缓 。 钴 头 处 所 受 轴 向 压力 与 钻 压 一 致 , 且 离 钴 头 

800 
700 
600 
500 


A1] 7J/kN 


T- 
> 
© 
er 
© 
= Wi 500 10001 p s 000 3 500 4 000 4 500 
局 (a) 轴 向 力 
co 
N 
e : 
300 转速 对 钻 柱 横向 振动 的 影响 
gar 是 不 同 转速 下 钻 柱 横向 变形 图 。 如 图 所 
© 200 ^r | -—- n | 
cC - | 
2 | 
O 


SE] 
| 


| [T 
90 r/min 


&K7 不 同 转速 下 钻 柱 横向 变形 图 


越 近 , 销 柱 所 受 压 力 越 大 。 在 图 6(b) 中 可 以 看 出 ， 
钻 柱 在 直 井 段 横向 变形 很 小 ,在 造 斜 段 与 直 井 段 交 
界 处 ,外 柱 受 井 壁 约束 ,在 该 位 置 发 生 较 大 变形 。 造 
斜 段 和 水 平 段 钻 柱 距 钻头 越 近 , 轴 向 压力 越 大 , 钻 柱 
变形 越 大 。 由 于 直 井 段 (井深 2700 m 以 上 ) 钻 柱 受 
拉 , 稳 定性 大 大 提高 ,其 横向 变形 只 有 0.1 mm。 
此 ,本 研究 主要 研究 造 斜 段 和 水 平 段 (井深 2 700 m 
以 下 ) 钼 柱 的 动力 学 特性 ,并 分 析 了 不 同 钻井 参数 下 
的 水 平井 钻 柱 非 线 性 横向 振动 行为 。 


(b) 横向 变形 


图 6 钻 柱 轴 向 力 和 横向 变形 图 


Fig.6 Axial force and lateral deformation 


示 , 不 同 转 速 下 钻 柱 横向 变形 的 形态 相近 ,这 说 明 近 
钻头 处 销 柱 横向 振动 受 转速 的 影响 较 小 。 


60 rmin 


3 goo ccn 
-3 


< 
78500 
Tani o 
Z5 6p $0 4/050 
o f 
S 72 E qaem 


120 r/min 


s 
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Fig.7 Lateral deformation of the drillstring under different rotation speeds 
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图 8 为 不 同 转速 下 钻 柱 节点 运动 状态 。 水 平 段 

钻 柱 受 重力 和 摩擦 力 的 影响 ,在 井 眼 右 下 角 旦 现 出 
圆 弧 摆动 的 运动 轨迹 。 由 于 近 钻 头 处 销 柱 受到 钻头 
影响 ,其 运动 轨迹 更 靠近 井 腿 中 心 ,所 以 距离 钻头 
越 近 ,节点 运动 轨迹 逐渐 呈现 出 对 称 的 趋势 。 同 一 
转速 下 ,由 于 近 销 头 处 横向 振动 更 大 ,水 平 段 钻 柱 
(节点 10, 距 钻头 30 m) 的 运动 状态 相 较 于 造 斜 段 
(节点 100 和 180, FE dl; 3k 300 m 和 480 m) 更 加 复 
杂 。 但 是 对 于 同一 节点 ,其 运动 轨迹 随 着 转速 的 改 
变 无 较 大 变化 ,可 见 转速 对 造 斜 段 和 水 平 段 钻 柱 横 
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向 振动 的 影响 很 小 。 已 有 研究 指出 ,在 无 弯曲 的 
纯 直 井 段 中 ,离心 力 随 着 转速 的 增加 而 逐渐 增 大 ,从 
而 增 大 了 接触 频率 、 切 向 摩擦 力 和 正 向 力 , 最 终 造 成 
钻 柱 横向 振动 加 大 。 然 而 ,由 于 钻 柱 偏心 距 小 
(2 mm) ,离心 力 远 小 于 接触 力 , 且 阻 尼 力 也 随 着 转 
速 的 增加 而 增 大 。 所 以 ,在 水 平井 钻 进 时 ,转速 对 外 
柱 横向 振动 的 影响 较 小 。 所 以 ,在 本 研究 钻 具 组 合 
下 , 当 转 速 为 30 ~ 120 r/min 时 ,改变 转速 对 水 平井 
钻 柱 横 向 振动 的 影响 较 小 。 
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不 同 转速 下 钻 柱 节点 的 运动 状态 


Fig.8 The motion state of the drillstring nodes under different rotation speeds 
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图 9 为 不 同 钻 压 下 钻 柱 横 向 变形 图 。 如 网 所 
AN ,在 本 研究 销 具 组 合 下 , 钻 压 为 5S ~60 kN 时 , 直 井 
段 销 柱 横向 振动 变化 很 小 ,但 造 斜 段 和 水 平 段 钻 柱 
横向 振动 有 较 大 变化 。 当 钻 压 低 于 40 kN BT, 钻 柱 
横向 振动 较 弱 , 横向 变形 次 数 较 少 , 当 钻 压 达 到 
60 kN 后 , 钻 柱 横向 振动 强度 明显 增强 , 钻 柱 频繁 产 
生 横向 变形 。 另 外 ， 由 于 近 钻 头 处 钻 柱 横向 变形 会 
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改变 销 头 侧 向 力 和 销 头 转 角 , 将 导致 钻头 运动 轨迹 
产生 偏差 ,增加 井 眼 轨迹 控制 难度 。 

图 10 为 钻 柱 节点 在 不 同 钻 压 下 运动 状态 。 可 
见 随 着 钻 压 增加 ,水 平 段 钻 柱 (节点 10) 的 运动 轨迹 
变 得 复杂 且 不 稳定 ,而 造 斜 段 钻 柱 (节点 100,180) 
的 运动 轨迹 变化 很 小 。 这 是 由 于 越 靠近 钻头 , 钻 柱 
所 受 压力 越 大 ,横向 振动 越 大 ,与 井 壁 接触 更 加 频 
繁 ,导致 钻 柱 稳定 性 急剧 降低 ,最 终 将 增加 近 钻 头 处 
钻 柱 失效 的 风险 。 
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到 9 不 同 钻 压 下 钻 柱 横向 变形 图 
Fig.9 Lateral deformation of the drillstring under different WOBs 
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图 10 不 同 钻 压 下 销 柱 节点 的 运动 状态 


Fig. 10 The motion state of the drillstring nodes under different WOBs 


3.3 扶正 器 对 钻 柱 横向 振动 的 影响 
图 11 为 不 同 扶正 器 安装 位 置 下 , 近 销 头 处 销 柱 
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的 横向 位 移 。 如 图 所 示 , ARR TRE B E CIR 
钻头 1.8、9、18、27 m) , 销 柱 产 生 很 小 的 横向 位 移 ， 


这 说 明 扶正 局 可 以 有 效 限 制 钻 柱 的 横向 位 移 。 在 该 
微 信 公 众 号 :应 用 力学 学 报 


第 3 Hj 


TE HABE P ,未 安装 扶正 器 时 , 销 柱 在 距 钻头 6m 处 
出 现 了 较 大 的 横向 位 移 ,而 扶正 天 一 般 只 能 限制 扶 
正 带 安装 位 置 附近 钻 柱 的 横向 位 移 ,所 以 当 扶 正 带 
安装 在 距 钻头 18 m 和 27 m 时 , 距 上 述 位 置 较 远 ,无 
法 有 效 限 制 近 钻头 处 销 柱 的 横向 位 移 , 其 作用 效果 
较 差 ; 当 扶 正 吉 安装 在 距 钻 头 9 m 处 时 , 距 上 述 位 置 
较 近 ,对 近 销 头 处 外 柱 横向 位 移 的 抑制 效果 更 好 。 
可 以 看 出 ,在 距 钻头 9 m 处 安装 扶正 器 不 仅 降低 了 
近 钻 头 处 销 柱 的 横向 位 移 , 且 钻 柱 横向 变形 次 数 远 
小 于 另外 4 种 情况 , 极 大 的 降低 了 近 钻 头 处 钻 柱 失 
效 的 风险 。 因 此 ,在 水 平 并 钻 进 时 ,为 了 防止 井 眼 轨 
迹 偏 移 , 减 小 近 钻 头 处 销 柱 的 弯曲 应 力 ,建议 安装 扶 
正 需 来 限制 钻 柱 的 横向 位 移 ,扶正 需 安 装 位 置 根据 


钻 柱 动力 学 特性 确定 。 
T= 000 
2 SERIO m 安 装 距 钻头 18m 安 装 
e 距 钻头 27 m 安 装 
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C5 图 11 不 同 扶正 器 安装 位 置 下 的 钻 柱 横向 位 移 


N Fig.11 The lateral displacement of the drillstring 
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-通过 上 述 研究 ,主要 结论 如 下 。 

〇 1 ) 在 水 平井 钼 进 过 程 中 ; 造 斜 段 和 水 平 段 钻 柱 
横向 振动 强度 较 大 , 且 越 靠近 钻头 , 钴 柱 振动 越剧 
烈 ,而 直 井 段 外 柱 横向 振动 较 小 ;水 平 段 和 造 斜 段 外 
柱 在 并 眼 右 下 角 呈 现 出 圆 弧 摆动 的 运动 轨迹 。 

2 ) 改 变 转速 对 水 平 并 全 井 段 铺 柱 横向 振动 影响 
很 小 ,改变 钻 压 对 直 并 段 钻 柱 横向 振动 的 影响 很 小 ， 
但 对 造 斜 段 和 水 平 段 钻 柱 横向 振动 有 较 大 影响 。 在 
本 研究 钻 具 组 合 下 , 当 钻 太 达 到 60 kN 后 , 钻 柱 横向 
振动 强度 明显 增强 。 因 此 水 平井 钻井 作业 中 ,在 满 
足 生 产 要 求 的 条 件 下 ,建议 优选 钻 压 来 降低 钻 柱 的 
横向 振动 ,提高 钻 柱 使 用 寿命 。 

3 ) 扶 正 器 可 以 使 水 平 段 钻 柱 在 井 眼 中 心 保持 稳 
定 , 且 扶正 器 安装 位 置 对 该 段 钻 柱 横向 位 移 的 影响 
较 大 。 所 以 ,在 水 平 并 销 进 过 程 中 ,建议 根据 钻 柱 动 
力学 特性 ,选择 合适 位 置 安装 扶正 器 来 限制 水 平 段 
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